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Isometries and Construction of Permutation Arrays 



Mathieu Bogaerts 
Abstract 

An (n, d) -permutation code is a subset C of Sym(n) such that the Hamming distance 
dn between any two distinct elements of C is at least equal to d. In this paper, we use 
the characterisation of the isometry group of the metric space (Sym(n), dn) in order to 
develop generating algorithms with rejection of isomorphic objects. To classify the (n, d)- 
permutation codes up to isometry, we construct invariants and study their efficiency. We give 
the numbers of non-isometric (4,3)- and (5,4)- permutation codes. Maximal and balanced 
(n, (i)-permutation codes are enumerated in a constructive way. 

I. Permutation arrays and permutation codes 

An (n, d) -permutation code of distance d, size s and degree n is a non-empty 
subset C of the symmetric group Sym(n) such that the Hamming distance between 
any two distinct elements of C is at least equal to d. The Hamming distance between 
two permutations 4>,ip E Sym(n) is defined as dni^,^) = \{i G {1, • • • , n} : <j){i) ^ 
vjj(i)}\. The sxn array A associated to a (n, d) -permutation code C = {(pi, . . . , <p s } of 
size s by Ay = 4>i(j) has the following properties: every symbol 1 to n occurs exactly 
in one cell of any row and any two rows disagree in at least d columns. Such an array 
is called a permutation array (PA) of distance d, size s and degree n. Permutation codes 
have first been proposed by Ian Blake in 1974 as error-correcting codes for powerline 
communication [2]. This application motivates the study of permutation codes having 
the largest possible size. Upper bounds for the maximal size fx(n, d) of a permutation 
code with fixed parameters n and d have been studied by many authors, see e.g. Deza 
[8], Cameron [6], and more intensively since Chu, Colbourn and Dukes [7],Tarnanen 
[18], and Han Vinck [1], [19]. 

A (n, d) -permutation code C of size s is maximal if C is not contained in a (n, d)- 
permutation code of larger size s' > s. Note that a (n, d) -permutation code reaching 
the maximal size /i(n, d) is necessarily maximal while the converse is not true. In this 
article, we will focus on the generation up to isometry of all maximal permutation 
codes with fixed parameters n and d. 

The paper is organised as follows: section [II] describes the isometry group of the 
metric space (Sym(n), dn), section |nl] contains improvements of methods to obtain all 
PA's of a given size. Section [IV] describes two algorithms for generating all maximal 
permutation codes up to isometry. 



II. Isometries 

A distance D on Sym(n) is called left-invariant (resp. right-invariant) if D((f),ip) = 
D(a<f>,aip) (resp. D^ip) = D(cj)a,ipa) ) for all a,(f),ip G Sym(n). A distance 
that is both left- and right-invariant is said to be bi-invariant. For any bi-invariant 
distance, the left multiplications l a : <p i— ► acj) and the right multiplications r a : cf) i— > 
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<j)a~ l are isometries. As noticed by Deza [9], any bi-invariant distance is invertible: 
D((f>,ip) = D^ -1 ,^ -1 ), or equivalently, the inversion 1, mapping each permutation 
onto its inverse, is an isometry. Let denote by 1Z (resp. C) the group of all right (resp. 
left-) multiplications and by X the group generated by the inversion 1. We will say 
that the distance D distinguishes the transpositions if there exists a constant c such 
that D(cf), ip) — c •<=>- (fti/j^ 1 is a transposition. 

In 1960, Farahat characterized the isometry group Iso(n) of the metric space 
(Sym(n), d H ) [11]. Since the Hamming distance is bi-invariant and distinguishes the 
transpositions, the following result generalizes the characterisation given by Farahat: 

Theorem 1 Let D be a bi-invariant distance distinguishing the transpositions on 
Sym(n) (n > 3), then the group Iso^ of isometries of (Sym(n), D) is (C x 1Z) x X, 
isomorphic to Sym(n) I 2. 

Proof: For any bi-invariant distance D distinguishing the transpositions, the isometry 
group of (Sym(n), D) contains the groups C, 1Z and X, generating the group (£ x 
1Z) x X. The pointwise stabiliser of the set containing the identity, the transpositions 
(1, k) with 2 < k < n and (1, 2, 3) in the isometry group Iso(n) is the trivial group, 
proving that there is no other isometries than (£ x 1Z) xi X. See [4] for developements. 

Every isometry t e Iso(n) can be uniquely written as ^r^i* with k — or 1, 
a, (3 G Sym(n). The action of a left multiplication l a on a given code corresponds to 
the permutation under a of the symbols appearing in the PA associated to the code, 
and the action of a right-multiplication r# is equivalent to the permutation under (3 of 
the columns of the PA. In other words, classifying permutation codes up to isometry 
is equivalent to classifying PAs up to permutation of their rows, their columns, their 
symbols and up to the inversion. 

III. Generate all (n, d) -permutation codes of a given size 

An efficient way to generate all permutation codes of a given size is described in 
[7] as a clique search. A graph is built whose vertices are the elements of Sym(n) and 
whose edges are the pairs of permutations at distance at least equal to d. A (maximal) 
(n, d)— permutation code is equivalent to a (maximal) clique in this graph. Since 
the isometry group Iso(n) is transitive on Sym(n), any (n, d) -permutation code is 
isometric to a permutation code containing the identity permutation Id. Assuming that 
the permutation code contains Id, the clique search is restricted to the neighbourhood 
of the vertex Id, that is, to the set of permutations that are at distance at least equal to d 
from the permutation Id. This improvement reduces the number of vertices to consider 
from n\ to Y^l=d wnere is the number of fixed-point-free permutations of 

k points . 

Using this method, we confirm the results of T. Kl0ve [13]: there are, up to isometry, 
exactly 7 (6, 5)— permutation codes of size 18. Note that the equivalence between PAs 
used in [13] doesn't take in consideration the inversion i. 



3 



The subgroups of Sym(n) can be used to construct large PA's. Starting form a 
subgroup G, Chu, Colbourn and Dukes used the orbits of L G =< lp : (3 E G > 
to simplify the construction of PA's using a clique search algorithm. More recently, 
Janiszczak and Staszewski [12] used the orbits a subgroup of /so(10) (more precisely, 
a subgroup XJ G of the conjugacies Uq =< Iprp : (3 E G >) in order to prove that 
^(10,9) > 49. Using other subgroups of Iso(n), we found other codes as e.g. an 
(7, 4)— permutation code of size 336 preserved by a subgroup of Iso(7) of order 7056 
[3]. 

IV. Generating all maximal permutation codes up to isometry 

In the past few years, techniques have been developped in order to construct PA's 
of large size. We adopt another approach, and for small values of n and d, succeed to 
classify all maximal permutation codes up to isometry. Note that for d = n, a maximal 
permutation code has size s = n and the PA associated to such a code is a latin square. 
The enumeration of latin squares up to permutation of the rows, columns and symbols 
has been widely studied and has been achieved for n < 11, see [16], [17]. As far as 
we know, the only classification of maximal (n, d) -permutation codes with d < n has 
been made by Torleiv Kl0ve for (6, 5)— permutation codes of size 18, as mentionned 
in section [TTTl 

In order to generate all maximal permutation codes with fixed parameters n and d, 
we used two isomorphism free generation algorithms, as introduced by McKay [14]. 
The first one is based on the backtrack principle, at every step the generated code 
is compared to the list of already encountered codes. We introduce here the notation 
Vd(C) for the set of permutations that are at distance at least equal to d of all elements 
in C. 



Algorithm 1 Generate a list L of all non-isometric codes 
Procedure GENBYLIST (C, L) 
if 3C" E L such that C ~ C then 

break 
end if 

LUC 
if C is maximal then 

return C 
end if 

for <f> E V d (C) do 

GENBYLIST( CU {(/>}, L) 
end for 



Initiated with an empty list and a code C consisting of a single permutation, algo- 
rithm \T\ produces exactly one representative of each permutation code up to isometry. 
It requires to compare permutation codes, motivating the developpement of invariants 
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that are efficient to test whether two permutation codes are isometric or not. 

We define the pair of quotient sets of a (n, d) -permutation code C as the ordered 

set (A(C),E(C)) where 

A(C) := {W,- 1 : hi/, e C} 

E(C) := {fV^^eC} 

Two pairs of quotients sets (A, E), (A', £') are equivalent if there exists a, (3 E 
Sym(n) such that aAa' 1 = A and /3E/3" 1 = E' or aAa^ 1 = E' and /3E/T 1 = A'. 

A famous code invariant is the distance enumerator polynomial of a code C, that 
enumerates the number of pairs of permutations wich are at distance k for < k < n. 
This invariant is weaker than the cycle index polynomial Qc( x ) of the permutation 

code C, defined using the p(n) conjugacy classes Ci, . . . , C p(n ) as Qc{x) = — 
where bj is the number of pairs of permutations ((f), 'ip) of C such that (pip' 1 E Cj. 
If the code C is a group, the cycle index polynomial of C coincides with the usual 
definition of the cycle index polynomial of a group. 

The set of occurences Occ(C) associated to a (n, d) -permutation code C is defined 
as Occ(C) := {ojj : 1 < i, j < n} where Oij is the number of permutations <p E C such 
that <J>(i) = j. If Occ(C) = {r}, then the permutation code is said to be r-balanced 
(see [10]). 

Given an (n, d) -permutation code C of size s, we define the graph (V(C), E(C)) as 
follows. The s + n 2 + 2n vertices are partionned into four sets V\ := : 1 < i < s}, 
V 2 := {cj : 1 < j < n}, V 3 := {sj : 1 < j < n} and V 4 := {dj : 1 < i,j < n}, 
(pi is the z-th permutation of C, Cj is the j-th column, Sj is the j-th symbol and Cjj 
is the (z, j)-cell of the permutation array associated to C. The edges of the graph are 
defined by the following rule: 0fe(i) = j if and only if the vertex E V4 is adjacent 
to the vertices <fi k E V\, Ci E V 2 and Sj E V 3 . 

Theorem 2 The pair of quotient sets, the cycle index polynomial and the set of 
occurences are invariants, but are not complete invariants. The graph (V(C), E(C)) 
of the code C is a complete invariant. 

The proof of this statement is routine check since the isometry group is described by 
Theorem Q] 

The graph (V(C), E(C)) can also be used to obtain a canonical form of the per- 
mutation code by a canonical labelling of the vertices of the graph (using e.g. routines 
of Nauty [15]). The canonical form automatically defines a canonical augmentation 
p(C) for any code C, and insures that the algorithm [2] generates all non-isomorphic 
maximal codes. 

Results of algorithm [T] are confirmed by the results of algorithm [2l there are, up to 
isometry, 61 (4, 3)— permutation codes (4 are maximal) and 9445 (5, 4)— permutation 
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Algorithm 2 Returns all non-isometric maximal codes 
Procedure CANONICAL-AUGMENTATION (C) 
if C is maximal then 

return C 
end if 

C^{CU{0}:0Gy d (C)} 
D^{Cn {*(£>) : £ G Stafe(C)} :DeC} 
for ZgD do 
choose K <E Z 

if 3i G Stab(K) s. t. t(C) = p(C) then 
CANONICAL-AUGMENTATION*^) 
end if 
end for 



TABLE I 

NON-ISOMETRIC MAXIMAL (5, 4) -PA'S 



Size 


Number of non-isometric (5,4)— PA's 


7 


1 


8 


25 


9 


36 


10 


46 


11 


18 


12 


10 


13 


1 


15 


1 


20 


1 



codes (139 are maximal). Table U details the numbers of non-isometric maximal PA's 
of distance 4 on Sym(5), of size 7 to 20. 

An r-balanced PA with n symbols has size exactly nr. We constructed all non- 
isometric r-balanced codes C(n,d) for small values of n and d. Table [n] shows the 
exact number M r (n,d) of non-isometric r— balanced (n, d)— permutation codes for 
small values of n, d and r. All our permutation codes and the code of our Magma [5] 
programs are available [3]. 

TABLE II 

Number of non-isometric t-balanced PA's 



r 


n 


d 


M r (n,d) 


2 


5 


4 


6 


2 


5 


3 


218 


2 


6 


5 


2799 


3 


5 


4 


1 


4 


5 


4 


1 
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